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Abstract 

The objective of this paper is to investigate how the breaking of manifest rotational invariance 
in light-front dynamics affects the binding energy and the form factors of the deuteron. To do 
this, we derive new light-front nucleon-nucleon one- and two-meson-exchange potentials, and use 
the potentials to solve for the deuteron wave function and binding energy. We find that includ¬ 
ing two-meson-exchange (TME) potentials partially repairs the broken rotational invariance of the 
one-meson-exchange (OME) potential. This is shown by a decrease in binding energy difference 
of different m states of the deuteron. We calculate the matrix elements of the electromagnetic 
current using the deuteron wave functions obtained from the OME and OME-I-TME potentials. 
Rotational invariance requires that the matrix elements satisfy an angular condition, but in light- 
front dynamics that condition is only partially satisfied. The current matrix elements from the 
OME calculation satisfy the angular condition better than the ones from the OME-I-TME calcu¬ 
lation. The matrix elements of the axial current satisfy the angular condition to the same extent 
regardless of which wave functions are used for the calculation. Einally, we hnd that at momentum 
transfers greater than about 2 GeV^, the breaking of rotational invariance causes less uncertainty 
in the computed deuteron form factors than do the uncertainties in the nucleon form factors. 

PACS numbers: 21.45.-tv, 03.65.Ge, 03.65.Pm, ll.lO.Ef 
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I. INTRODUCTION 


Recent experiments at Thomas Jefferson National Accelerator Facility have measnred 
the A{Q‘^) strnctnre fnnction of the denteron for momentnm transfers np to 6 (GeV/c)^ [|l|, 
and measnrements for B{Q‘^) are planned. Eventnally, even higher momentum transfers will 
be achieved. At such large momentum transfers, a relativistic description of the denteron 
is required. Even at lower momentum transfers, a relativistic description is important to 
understand the details of the form factors. In addition, incorporating relativity is important 
for the denteron wave function to transform correctly under boosts to large momentum, 
which is important for calculating form factors. 

One approach that gives a relativistic description of the denteron is light-front dynamics. 
The subject of this work is to investigate the consequences of combining light front dynam¬ 
ics with various nuclear models to calculate bound state wave functions, and use them to 
calculate the denteron form factors. 

The utility of the light-front dynamics was first discussed by Dirac p|. Light-front dy¬ 
namics makes use of the light-front coordinate system, where a four-vector is expressed 
as = {x~^ , x~ , x^ , x'^) , with x^ = ± This is simply a change of variables, but an 

especially convenient one. Using this coordinate system and defining the commutation rela¬ 
tions at equal light-front time {x~^ = ^lf), we obtain a light-front Hamiltonian |, |^. We 
use Hamiltonian in the light-front Schrodinger equation to solve for bound states. There are 
many desirable features of the light-front dynamics and the use of light-front coordinates. 

First of all, high-energy experiments are naturally described using light-front coordinates. 
The wave front of a beam of high-energy particles traveling in the (negative) three-direction 
is defined by a surface where x~^ is (approximately) constant. Such a beam can probe the 
wave function of a target described in terms of light-front variables BBH: the Bjorken x 
variable used to describe high-energy experiments is simply the ratio of the plus momentum 
of the struck constituent particle to the total plus momentum (p+) of the bound state. 

Secondly, the vacuum for a theory with massive particles can be very simple on the light 
front. This is because all massive particles and anti-particles have positive plus momentum, 
and the total plus momentum is a conserved quantity. Thus, the naive vacuum (with p'^ = 0) 
is empty, and diagrams that couple to this vacuum are zero. This greatly reduces the number 
of non-trivial light-front time-ordered diagrams. 

Thirdly, the generators of boosts in the one, two, and plus directions are kinematic, 
meaning they are independent of the interaction. Thus, even when the Hamiltonian is trun¬ 
cated, the wave functions will transform correctly under boosts. Thus, light-front dynamics 
is useful for describing form factors at high momentum transfers. 

Finally, it is easy to perform relativistic calculations using light-front dynamics. This is 
partly due to the simplicity of the vacuum, and partly due the the fact that, with light- 
front variables, center-of-momentum variables can be cleanly separated from the relative 
momentum variables. This allows us to write relativistic equations which have the simple 
form of a non-relativistic Schrodinger equation. 

One serious drawback of light-front dynamics is that rotational invariance is not manifest 
in any light-front Hamiltonian [^. This is a result of selecting a particular direction in space 
for the orientation of the light-front. 

An untruncated light-front Hamiltonian will commute with the total relative angular 
momentum operator, since the total momentum commutes with the relative momentum. 
Thus, eigenstates of the full Hamiltonian will also be eigenstates of the angular momentum. 
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However, as mentioned earlier, a Fock-space truncation of the light-front Hamiltonian results 
in the momentum operator four-vector losing covariance under rotations. Hence and the 
truncated Hamiltonian do not commute and this implies that the eigenstates of the truncated 
Hamiltonian will not be eigenstates of the angular momentum. 

How will this violation of rotational invariance affect physical observables? One way to 
observe this violation is to note that on the light front, rotational invariance about the z-axis 
is maintained. This allows us to classify states as eigenstates of J 3 with eigenvalues m. We 
compare the binding energy of deuteron states (which have j = 1 ) with different m values. 
If the Hamiltonian were rotationally invariant, the energies should be the same; the breaking 
of rotational invariance causes the energies to be different [^. 

Another symptom of the breaking of rotational invariance is that the angular condition 
(a relation between the matrix elements of the current, dehned in section [IV A 4|) for the 
deuteron current is not exactly satished 0 0 0 , [I^ , [T^ . This means that different 


prescriptions for calculating the deuteron form factors from the deuteron current will in 
principle give different results when light-front dynamics is used. This dependence on the 
prescription used has caused concern about the validity of applying light-front dynamics to 
calculate form factors. We address both the breaking of binding-energy degeneracy and the 
angular condition in this paper. 

One notable feature of this calculation is that it is done entirely with light-front dynamics. 
The covariant Lagrangian generates light-front potentials, which generate light-front wave 
functions, which are used in a light-front calculation of the deuteron current and form factors. 
This is different from other approaches which use deuteron wave functions calculated from 
equal-time dynamics, then transformed to the light front []^ [^, |^, |K, 




We start by introducing a model Lagrangian for nuclear physics which includes chiral 
symmetry in section H. The methods introduced in Refs. |]T 8 | , p!9| are generalized for use 
with this nuclear model. The Hamiltonian is derived and used to calculate new light-front 
nucleon-nucleon one-meson-exchange (OME) and two-meson-exchange (TME) potentials. 
The notation and conventions dehned in the appendix are used extensively in this section. 
We have some freedom in how to choose the TME potentials, and we consider several 
different choices. 

In section m the OME and TME potentials are used to calculate the deuteron wave 
functions and binding energies. We hnd that by including the TME potentials for the 
calculation of the deuteron, the binding energy degeneracy is broken by a smaller amount. 

The wave functions obtained in section m are used in section E to calculate the electro¬ 
magnetic and axial form factors of the deuteron. Although rotational invariance demands 
that there be only three independent components of the deuteron current, the light-front 
calculation of the deuteron current results in four independent components. This is a result 
of the lack of manifest rotational invariance on the light front. There are several prescrip¬ 
tions for choosing which deuteron current component should be eliminated, and in principle 
this choice will affect the form factors. We attempt to hnd currents that transform correctly 
(or well enough) under rotations so that the choice of “bad” component does not matter too 
much. 

We discuss the results of that search in section |V[ We hnd that for most of the currents, 
the angular condition does not depend strongly on which potential is used to calculate the 
wave function. The only exception to this is that part of the electromagnetic current which is 
multiplied by the isoscalar Fi nucleon form factor satishes the angular condition much better 
when using the wave function calculated with the OME potential than with wave function 
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with other potentials. We also find that the major uncertainty in the calculated deuteron 
form factors at momentum transfers greater than 2 GeV^ is due not to the prescription used 
to determine the form factors from the current, but instead is from the uncertainties in the 
nucleon form factors. 


II. REALISTIC NUCLEAR MODEL 


We consider a nuclear Lagrangian which uses (tt, a, p, uj, rj, and h) for the nucleon-nucleon 
interaction. This is used to calculate a new light-front nucleon-nucleon potential (LFNN). 
This model is an extension of the light-front model used by Miller and Machleidt |^. A 
new feature of this model is that light-front energy dependent denominators are used in the 
potentials; the denominators used in Ref. are energy independent. 

Note that most of the material covered here can be considered an extension of the work 
done with the Wick-Cutkosky model in the Refs |^. This allows us to build upon the 
previous results. 


Model and Formalism 


Our starting point is a nuclear Lagrangian 0 which incorporates a non-linear chiral model 
for the pions. The Lagrangian is based on the linear representations of chiral symmetry used 
by Gursey [^. It is invariant (in the limit where mjr ^ 0) under chiral transformations. 

The model prescribes the use of nucleons '0 (or V'O ^rid six mesons: the vr, 6 (also known 
as the ao(980)), a (also known as the /o(400 — 1200)), p, p, and oj mesons. The coupling of 
each meson to the nucleon is governed by the combination of the meson’s spin and isospin. 
The TT and p are pseudoscalars, the p and uj are vectors, and the 5 and a are scalars. Under 
isospin transformations, the tt, p, and 5 are isovector particles while the p, uj, and a are 
isoscalar particles. 

The use of scalar mesons is meant as a simple representation of part of the two-pion- 
exchange potential which causes much of the medium range attraction between nucleons 
It can also be interpreted as the effect of fundamental scalar mesons 

The Lagrangian C is based on the one used in Refs. [|, |2^, |2^ 
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It is given by 


C = 




mi 1 

-^p>^ ■ p^ - + -^UJ^^UJ^ 


+^fTT{d,U + \mlfTT{U + Ut - 2) 

+^id^ad^a - m^a^) ■ d^5 - + ]^{d^pd^p - m^p^) 

+'ip' - ppp^ ■ T - g^ujp) -U{M + g^a + gs5 ■ T + ig^'y^p) iIj', (1) 


where the bare masses of the nucleon and the mesons are given by M and ruo where a = 
n,p,a,6, p,u!. We have defined = d^V'" — d'^V^ for V = p,uj. The unitary matrix U 
can be chosen to have one of the three forms Uj. 


Ui = U2 = 


l + ip 5 T ■ 7r/2/ 

1 - iysT • 77/2/’ 


f /3 = Ul - 7r7/2 -h i'y^T • tt//, 


( 2 ) 
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which correspond to different definitions of the fields. Note that each of these definitions 
can be expanded to give 


U 


1 + *75 


T ■ TT 

“T 




( 3 ) 


In this work, we consider at most two meson exchange potentials, so we consider U to be 
defined by Eq. (||). 

In the limit where > 0, this Lagrangian, is invariant under the chiral transformation 




„i 75 X-a 




U — > g-^^ST-a JJ 


( 4 ) 


In this model the other mesons are not affected by the transformation because they are not 
chiral partners of the pion. This is in contrast to the Lagrangian given in Refs. @.0. where 
the mass and scalar interaction terms for the nucleon were written as MIJ + ggcj) instead of 


B. Non-interacting Nncleon-Nncleon Theory 


The light-front Hamiltonian is derived from this Lagrangian using the same approach 
used in Refs. |T^, |T9[ , the approach used by Miller 0 and others ^ . The basic 

idea is to write the light-front Hamiltonian {P~) as the sum of a free, non-interacting part 
and a term containing the interactions. We consider the free part first. 


1. Free Field Expansions 


The solutions for the free fields are similar to those obtained by using equal-time dy¬ 
namics. In fact, the solutions are formally related by a change of variable, and so the most 
obvious difference between the two is due to the Jacobian. The field equations have the 
general form (when Lorentz, spinor, and isospin indices are suppressed) of 


a{x) = 


d^k±dk+e{k^ 




(27r)3/2y^ 

where a = tt, r^, a, 6, p, to, Note that in the exponentials. 


k^x^ = - (^k~^x + k x~^^ — k± 


■ Xx. 


In particular, the solutions for all the mesons and the nucleon field are 
d‘^k±dk^6{P 


TZ{X} = 

rj{x) = 

5(x) = 


(27r)3/2^/^ 

d?kxdk+e{k-^ 


{2'KfPV^ 

d?kxdk+e{k+) 




( 5 ) 


( 6 ) 

( 7 ) 

( 8 ) 
( 9 ) 
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a{x) = 

J (27r)3/V2A;+ ^ ^ ^ J ’ 

(10) 

p^{x) = 

J (27r)3/V2A;+ ^ ^ ^ 

(11) 

OJ^{x) = 

/ i: s) [aUk, + aiik, , 

(12) 

%Ij{x) = ■ 

J (27r)3/V2A;+ 



X E E [u{k,X)b{k)e-^’‘''^>^+v{k,X)d\k)e+^^"^-]xt,. 

(13) 


A=+,— ts— 


The polarization vectors are the usual ones. The most general of the commutation relations 

is 

aa,i{k,s),alj{k',s') = 6a,pSij6s,s'S^'^’^\k - k'), (14) 

where a, i, and s denote the meson type, isospin, and spin. The anti-commutation relations 
are 

{b{k,X),b\k',X')} = [d{k,X),d\k',X')} = 6x^x'S^^’^\k-k'). (15) 

All other (anti-)commutation relations vanish. The spinors are normalized so that 
u{p, X')u{p, X) = 6x'x- For more information on the dehnition of the spinors, see the Ap¬ 
pendix. 

2. Non-interacting Hamiltonians 

The general form of the non-interacting Hamiltonian for each meson is 

/ 2 2 

Pq“(q;) = f d‘^k±dk~^9{k^)al,{k)aa{k) ■ (16) 

J K' 

For the vector mesons (p and u), there is an implicit sum over the meson spins. Explicitly, 
this means that for vector mesons a|/(fc)av'(fc) — s)av{k, s). Likewise, the 
sum over the isospin of the isovector mesons (tt, <5, and p) is implicit. The sum over isospin 
can be made explicit by writing a\{k)ai{k) Z]i=i, 3 j(^)o/,i(fc)- The non-interacting 
Hamiltonian for the nucleons has a similar form as well. 



These equations are what one expects, since a free particle with momenta k± and k~^ has 
light-front energy k~ = —^ 7 :+—. 
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C. Interacting Nncleon-Nncleon Theory 


The interaction Hamiltonians are derived from the Lagrangian in Eq. (|l]) using the tech¬ 
niques presented in Refs. O . However, there are some additional complications due to 


the structure of the interactions. 

One complication is that the chiral coupling of the pion field to the nucleons through 
the U matrix generates vertices with any number of pions. This is addressed simply by 
expanding the U matrix in powers of j, and considering the interaction Hamiltonians order 
by order. 

Another complication is due to the fact that both the vector mesons and the fermions 
have components which depend on other components of the field [^, . Vector 

meson fields have four components, but only three degrees of freedom. Likewise, fermion 
fields have four spinor components, but only two degrees of freedom. When the dependent 
components are expressed explicitly in terms of the independent components, we obtain 
new (effective) interaction Hamiltonians for instantaneous vector mesons and fermions. A 
complete derivation is given by Miller in Ref. p and earlier workers cited therein. We 
illustrate only the main points of the derivation here. 


1. Expanding the Pion Interaction 


We start by Taylor-expanding U in powers of after which the derivation of the 


one- 


meson-interaction Hamiltonian P'j ^ is straightforward. (The prime indicates that it is in 
terms of "0', not 'ip. We derive the expressions for PfR in the section |H C 2| .) The result is 


P'i,i = j R (x) + 95di{x)Ti + ga(j{x) 

+ 9Ahb)TiMx) + 9xihb)x{x) 


(18) 


We have dehned a dimensionless coupling constant = y. To save space and to generalize, 
we dehne 


if a is a pseudoscalar meson (tt, rj) 
1 if a is a scalar meson [6, a) 

7 ^^ if a is a vector meson (p, a;) 

Ti if a is an isovector meson (tt, S, p) 

1 if a is an isoscalar meson (p, a, u) 


and denote the meson fields by This allows us to write 


Ppi = H / d‘^x±dx ^jJ\x)gaTc,Ta^a{x)^jJ'{x), 

a=7r,7],a,S,p,uj 


(19) 

( 20 ) 


( 21 ) 


where the appropriate sums over the meson indices are implicit. 

The next step is to derive the two-meson-interaction Hamiltonian which arises from chiral 
symmetry, 


P 


1,2c 


= / d?x±dx Ip (x) 


2M 


TiTj7li{X)Tlj{X) 
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M 

9tt9x 

~w 


TiTTi{x)x{x) 'll)\x) 


= E 


(y.=7r,7],a,S 


9 - 1^90 

~w 


Sa / (fxxdx i!{x)\r„T„^Ax)\\raTa^^{x)]i^'{x), 


( 22 ) 

(23) 


where Sa is a symmetry factor, equal to | when a = vr, and 1 otherwise. When the 
contact interaction is used to calculate diagrams, an additional factor is picked up for the 
TTTT contact term (due to indistinguishability) which cancels the symmetry factor Note 
that this contact interaction involves only scalar and pseudoscalar mesons. 


2. Elimination of Dependent Fermion Components 

We are now ready to express the dependent components of if' in terms of the independent 
components, and address the problem of instantaneous nucleons. The generation of instan¬ 
taneous interactions is a general feature of theories with interacting fermions in light-front 
dynamics. This allows us to use a simplihed model to demonstrate how these instantaneous 
nucleons arise. In particular, we want to postpone the discussion of the complication intro¬ 
duced by the vector mesons until the next section. To this end, we choose to remove all 
mesons except the a from the Lagrangian given in Eq. (|l]). (The a is chosen since it has the 
simplest coupling to the nucleon.) The equation of motion for the nucleons is then 

ipfj' = (M + (24) 

Applying the projection operators A± = (dehned in the Appendix) to Eq. (|^) splits 

it into two equations, 

= [cx± ■ p± + P{M + g^a)]^jJ'_, (25) 

= [cx± ■ + g„a)]'il:'^, (26) 

where V'i = ^±4’'■ This split is useful because in Eq. (|^), all four components of the 
nucleon held are interrelated, while in Eqs. (|^) and (|^) , the two components of '0+ are 
related to the two components of ip'_, and vice versa. 

First, notice that Eq. (^) involves d~, a dynamic operator in light-front dynamics. 
Dynamic operators should be avoided since they involve the interaction, and are therefore 
complicated. We use Eq. (|^) to avoid that complication and relate the components of 
xjj. Secondly, to keep the relation as simple as possible, we do not attempt to invert the 
spinor matrix on the right side of Eq. (pS]). Requiring that the equation for the dependent 
components be both a kinematic equation and simple equation forces us to choose as 
the independent components. The dependent components, ipf, are dehned by 

i’- = ^ [q:±-P n+ /?(M-F5(^a)]'0+. (27) 

Notice that the dependent components consist of a non-interacting part and an interacting 
part. We separate these parts by dehning -0 (without a prime) to be the free nucleon held, 
and ^ to be the part of xf' that is due to interactions. So 

if' = V' + C, 
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(28) 






where 


^ = = ■\p{9ct(t)'iP+ = ^{g„a)'ip+. (29) 

2p+ 

This allows us to write 

^' = 'fp + ^{9a(T)^p+ (30) 

= + ^i9a(^)pj' (31) 

The last equation follows is obtained from noting that (q"*")^ = 0, which implies that = 

7+'0+. 

Plugging Eq. (^) into Eq. (|T 8 |) and removing all mesons except the a meson, we obtain 


-P 1,1 — ^ 1,1 + -P/,2) 

Pfi = / cfx±dx~ip{x)g^a{x)ijj{x), 


(32) 

(33) 


P /2 = / d'^x^dx ijj{x) 


'1 

9^(^{x)—g„a{x) 


pj^x) 


(34) 


We interpret the factor as a type of nucleon propagator that joins any two meson 
interactions (although having two of these propagators adjacent to each other causes the 
interaction to vanish since ( 7 +)^). Because this propagator does not allow for an energy 
denominator (as it is already between two potentials), ^ is called an instantaneous prop¬ 
agator. 

Thus, when constructing the diagrams for the light-front potentials, we must also include 
instantaneous propagators for the nucleons in addition to the usual propagators. 


3. Elimination of Dependent Vector Meson Components 


Like the nucleons, the vector mesons have a dependent component that contains interac¬ 
tions and must be eliminated. This process is complicated somewhat by the fact that the 
dependent nucleon components must be eliminated at the same time. The salient points 
of the combined elimination of the dependent nucleon and vector meson components are 
discussed in detail by Miller 0. 

The result is that the vector meson held must be redehned and an instantaneous vector 
meson propagator is generated in addition to an instantaneous nucleon propagator. However, 
when the nucleon-nucleon potential is calculated, the redehnition of the vector meson held 
exactly cancels the contribution of the instantaneous vector meson. The result is that the 
potentials can formally be calculated using the original vector meson held. 

In this work, we use that result to simplify our derivations of nucleon-nucleon potentials 
by formally using the nai'ive form of the vector meson held. Thus, we hnd an interaction 
Hamiltonian similar to the one shown in Eq. (^), 



Ppi + P, 


/, 2 ) 


(35) 
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( 36 ) 


P/,1 ^ j <Px_i_dx i\x)(joTJTa.^^{x)i^{x), 

a=n,r},a,5,p,uj 

Pi~2 = J d'^x±dx~'ilj{x)[ga^Ta^T„^<S>o.^{x)]^ 

c^l 5Ck2=7r, 77 ,(J,(5,p,cj ^ 

[ga 2 Ta 2 Ta 2 ^a 2 {x)]^^^{x). (37) 


We may continue to interpret as an instantaneous nucleon propagator, since in the 
derivation in Ref. ^ it is clear that the potential vanishes when there are two adjacent 
instantaneous propagators. In Refs. |^, sign of the coupling of the vector mesons in the 
equations equivalent to Eq. (^) has the wrong sign; the coupling of the mesons in Eq. (| 
must be the same as in Eq. 


Also note that in principle the same prescription has to be applied to the contact interac¬ 
tion, although to the order of two mesons, this simply has the effect of removing the primes. 
We hnd that 


P; 


1,2c 


= H J d'^x^dx ‘ip{x)[g^r^T^^^{x)][gaT^T^^^{x)]'ijj{x). (38) 

a=7r,77,cr,(5 


Note that Pj 2 c Pi 2 have forms that are very similar. In fact, we can obtain Pj 2 c from 
Pf 2 by making the following changes: 

1. Replace ^ with 

2. Replace cxi with tt. 

3. Restrict the sum on 02 to vr, g, a, S. 


4 . Interaction Hamiltonians in Momentum Space 


We look at the matrix element of the interaction Hamiltonians between initial and hnal 
states given by |fcj,Aj,rj) = 6l(/cj, Aj)^^-^ |0) and {kf,Xf,Tf \ = (0|6(fc/, A/)^!^ where both the 

bra and the ket have units of We hnd that 


PfAfn) 


gc 2 M r d'^q_Ldq^e{q+) 
a=n^,s,P,^ 2(271)^kpt ' ^2{2nf^ 

u{kf, Xf)Tau{ki, Xi)xl^TaXn 


d?Xi_dx 




.(q)e" 


Hai{q)e 


+iq'^Xfj, 


(39) 


where q = {q~^,q±) and there are implicit 6*-functions on for each particle to ensure 
that the light-front energy is positive. These are occasionally suppressed to simplify the 
equations. We include a meson-nucleon form factor to phenomenologically account for 
the fact that the mesons and nucleons are composite objects. 

Now, we evaluate P~ a.t x^ = 0 and use 

J = 2{27r)^6^^'P{kf - h - q), (40) 
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to write 




a=TT,ri,a,5,p,Lij \J 2(27r)^ (J~^ 

aa{q) 0 {kj -kt) + ai{q)e{-{kj - kf))] F^{q). 


(41) 


We define q = sign(A;^ — kf){kf — ki), which ensures that the q~^ of the meson is positive. 

Next we consider the two-meson interactions. First take the interaction with the instan¬ 
taneous propagator given by Eq. (|37|). Plugging in the held dehnitions gives 


= E 


f d‘^qi±dqt6{qt) f d‘^q2±dq2 0{q^ 


ai,a2=.,,,.,5,p,.2(27r)3^fc^^ J (27r)3/2^ 

\f)To,^'j+To,^u{ki, Xi)xlfTa^To,,XTi 




I + ai^{q 2 )e-^^P^^'^\ F^,{q 2 ) 

+ at^(qi)e+<"''] F^Mi) 

Y, A;)r^ 27 +I^( fc^ 

ai,a2=n,r},(T,S,p,u} 2{2 tt)^ k~j kf 


xlfT^2TaiXn 


d^km±dk+6{k+) 


tta^ikf - km)0{kj - k+) + ai^{km - kf)e{k+ - kj)] Fo^M-i) 
O^aiikm ki)0(^k^ ^m) Pciii.Ql') 


(42) 


(43) 


Note that the momenta qi and q 2 are implicitly functions of the momenta kf, km, and fcj. 

We could also write out the contact interaction given by Eq. ( pH]) in momentum space, 
but it is related to Eq. (|^) by replacing with ^ and restricting the allowed values 

of the a’s. 


D. Feynman Rules for Nucleon-Nucleon Potentials 


Now that we have the one-meson- and two-meson-exchange expressed in momentum 
space, we are now ready to write out the Feynman rules for diagrams in our model. For 
simplicity, the only diagrams considered are those where a meson emitted by one nucleon 
is absorbed by the other nucleon. Since we are only interested in the two-nucleon to two- 
nucleon potentials, we follow the same approach as outlined in Ref. to derive the rules. 
We denote a “normal” nucleon propagator by a solid line with an arrow, an instantaneous 
nucleon propagator by a solid line with an stroke across it, mesons of all type by a dashed 
line, and energy denominator terms by a vertical, light, dotted line. 

1 . Overall factor of ^ _ 


2(27r)= 


iU+ P.+ ju~t" p."t" 


2. Usual light-front rules for p± and momentum conservation. 


11 

































3. Factor of for each internal line, including any instantaneous nucleon lines. 


4. A factor of 




for each energy denominator. 


5. Each meson connects the two nucleons, and each end of the meson line has a factor 
of gaTaTaFa{q). The indices of the isospin factors on each end of the meson are 
contracted together. The Lorentz indices of the gamma matrices are contracted with 
—g^'^ for the vector mesons. 


6 . For each contact vertex, multiply by a factor of -T. If the vertex is a vr — vr vertex, 
symmetrize the = TiTj by replacing it with Sij. 

7. Factor of u{k, X)u{k, A) for each propagating nucleon and ^ for an instan¬ 

taneous nucleon. 


8 . Integrate with 4M^ / internal momentum loops. 

9. Put the spinor factors for nucleon 1 and 2 between uFs and the isospin factors between 
the initial and hnal state isospin. 


From this list, it is useful to summarize what needs to be done to convert a graph with 
an instantaneous nucleon to one with a contact interaction: 

1. Replace with 

2. If both mesons are pions, replace T^Tj with 6 ij. 

These rules make it easy to write down what various potentials are. 


E. Nucleon-Nucleon Potentials 


The meson exchange potentials have the same basic form as in Refs. [18 


|], however we 

must include the contact interaction and instantaneous nucleon propagators for the nuclear 
model used here. First, we discuss how to include the contact diagrams from the standpoint 
of the Bethe-Salpeter equation, then we begin to calculate the light-front potentials. 


1. The Bethe-Salpeter Equation and Chiral Symmetry 


The kernel of the Bethe-Salpeter equation 
richer than the one presented in Refs 


40 


for this nuclear model is 
I9[| for the Wick-Cutkosky model. This is due mainly 


to the presence of the contact interactions which are generated by the chirally invariant 
coupling of the pion to the nucleon. Several of the lowest-order pieces of the full kernel K 
are shown in Fig. |l|. (Note that for Feynman diagrams, it is useful to combine the “normal” 
nucleon propagator with the instantaneous nucleon propagator |]4l| , |4^ , |43[| , and denote the 
combination with a solid line.) 

As discussed in Refs. ||I8|, p!9|, each of these Feynman diagrams is covariant. This means 


that we may choose any of the diagrams from K to construct a new kernel K\ and the 
inhnite series of potential diagrams physically equivalent to K' will also be covariant. 
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In practice, this means that when deciding which two-meson-exchange potentials to in¬ 
clude for calculating the deuteron wave function, we may neglect the crossed diagrams. 
Although including only the box and contact two-meson-exchange diagrams may affect the 
exact binding energy calculated, we should hnd a partial restoration of rotational invariance. 
We reiterate that the focus of this work is on understanding the effects of the breaking of ro¬ 
tational invariance and how to restore it; our goal is not precise agreement with experimental 
results. 

We also want to keep the potentials chirally symmetric as well. Whereas Lorentz sym¬ 
metry is maintained by using a kernel with any Feynman diagrams (with potentially arbi¬ 
trary coefficients), chiral symmetry relates the strength of the tttt contact interaction to the 
strength of the pion-nucleon coupling. 

Chiral symmetry tells us that for pion-nucleon scattering at threshold, the sum of the 
time-ordered graphs approximately cancels [^]. Furthermore, upon closer examination, we 
hnd that all the light-front time-ordered graphs for the scattering amplitude vanish except 
for the two graphs with instantaneous nucleons and the contact graph. These graphs are 
shown in Fig. ^ Using the Feynman rules, and denoting the nucleon momentum by k and 
the pion momentum by g, we hnd that 


Mu = C 


TiTj 


Mx = C 


2 {k+ + q+) 

TjTi 


.-4 




u{k)Mu{k), 

u{k')'y~^u{k), 


Me = C^u{k')u{k), 

where the factors common to all amplitudes are denoted by C. 
For threshold scattering, we take k = k' = M and q = q' = 


Mu = C 
Mx = C 
Me = O 


/ ^i,j T 
2{M + m^) 

/ i^i,j,kTk 

2(M — m-n) 

/ ~^i,j 

M ’ 


(44) 

(45) 

(46) 

In that limit, we hnd 

(47) 

(48) 

(49) 


where C = Cu{k')u{k). In the limit that tUtt —>■ 0, the sum of these three terms vanishes. 
The term in these equations proportional to is the famous Weinberg-Tomazowa term 

mm. 


The fact that the amplitudes cancel only when the contact interaction is included demon¬ 
strates that chiral symmetry can have a signiheant effect on calculations. In terms of two- 
pion-exchange potentials, this result means that the contact potentials cancel strongly with 
both the iterated box potentials and the crossed potentials. This serves to reduce the 
strength of the total two-pion-exchange potential, which should lead to more stable results. 

However, since we do not use the crossed graphs for the nucleon-nucleon potential, we 
must come up with a prescription which divided the contact interactions into two parts 
which cancel the box and crossed diagrams separately. We do this by formally dehning two 
new contact interactions, so that 


Meu 


M 

2{M + rriT,) 


Me, 


(50) 
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Mc^ = Mc-Mc„. 


( 61 ) 


With these dehnitions, we hnd that at threshold and in the chiral limit, 

Mu + Mc^ = Q, (52) 

Mx + Mc^ = 0. (53) 

This indicates that we can incorporate approximate chiral symmetry withont inclnding 
crossed graphs simply by weighting each graph with a contact interaction by a factor of 

M 

2 (Af+m 7 r) * 


2. OME Potential 


The one meson exchange (OME) potential connects an initial state with two nncleons to 
a hnal state with two nncleons, has one meson in the intermediate state, and has the meson 
emitted and absorbed by different nncleons. With these restrictions, along with the fact 
that in light-front dynamics the interaction does not allow for particles to be created from 
the vacnnm, we hnd that there each meson has only two diagrams for the OME potential. 
These diagrams are shown in Fig. ^ 

The Feynman rnles derived in the previons section are used to derive the potential for 
these diagrams. We factor out an overall factor of common to all 

the two-nucleon potentials and suppress it from now on. Then we get 


Tome, a 


“ (27r)3 ’ Al0w(fc2/, A2/)r„n(fc2„ \2r)FM 

■_0 {kt} - kti) _ 

Xktf - kti){P- - ku - k^f) -ml- {ku^x - kif^±y 
_ Ojkti - ktf) _■ 

{kti - ktf){P~ - k^f - k^i) - ml - {k^^x - fci/,±)2_ ‘ 


(54) 


In the scattering regime (P = -|- /c 2 / = Eq. (|5^) agrees with Eq. (4.7) in 

Ref. i, after taking into account the difference in spinor normalization (we use uu = 1). 
To simplify the potential, we dehne 


a = 


0{k+f - 4+) X (4+ - 4+)(F- - if, - ky) + 


mi - kt,) X («:+ - kt,)(P- - K, - 

h = 2ki^±kf^j_, 




so 




0ME,Q; 


da'^a,! ■ 


fM 


Aij)rQ,'u(fcij, \\i)u{k2f 1 ^2i) 


( 55 ) 

(56) 


(57) 


(27r)3 [a-ml + bcQs[<j>f - 4>i)] 

Now we consider the precise form to use for the meson-nucleon form factor. We assume 
that it has a n-pole form ||^, so that the denominator of the meson-nucleon form factor has 
the same form as the denominator of the potential in the scattering regime. In particular, 
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playing the role of rria for the form factor. For simplicity, we declare that the denominator 
of the form factor always has the same form of the denominator of the potential. Thus, 


Fa{q) 


FI 


m: 


- A2 + hcos{(j)f - (j)i) ^ 


(58) 


Inserting the explicit expression of for the meson-nucleon form factor into the potential 
results in 




0ME,Q; 


/a-*-O', 


1 ■ T’a ,2 


( 2 vr) 


(A^ - <) 


2 \ 2nn 


[a - ml + bcos{<ijf - 0,)] [a - + ;>cos(.^/ - 0;)]^"° 


(59) 


Now, since the light-front potentials have exact rotational invariance about the z-axis, 
they conserve the Jz quantum number m. Thus, these potentials connect only states with 
the same value of m. This means that, in general, light-front potentials may be written as 


OO 

m=—OO 


(60) 


where is the potential in the magnetic quantum number basis. This relation can be 
inverted to obtain in terms of (pi), 


vm 

'^OME,a 


9lTa,i-T, 


0,2 / A 2 


(A^ - 


2 \ 2na 


(27r)3 

d(pf i^^^ ['^(^1/) ^lf)^aU{kii, Xli)u{k2f, X2f)TaU{k2i, A2i)]^_Q 


27r 


[a - m2 &cos(0/)] [a - A2 &cos(0/)] 


2?t,q 


(61) 


To continue on, we need to get an expression for (pf dependence the uu matrix elements. 
These are calculated in Ref. |^. Summarizing, we can write 

(62) 


“(fci/,Ai/)r„ti(fc„,Ai.)ti(fc2;,A2/)r.,M(fc2i, Aji) = y)c„(r„,r„)e‘'*', 


where the C depends implicitly on all the variables on the left-hand side except (pf and (pi. 
Thus, 

t _ ™ 2 '\ 2 nQ, 


d(pf 


(27r)3 


27r [a — m^ -I- 6cos(0/)] [a — A^-1-&cos(0/)]^”'“ ^ ^ 

Substituting / for the 


The cosine integral is denoted by I, and is calculated in Ref. 
integral, we obtain 

,2 A2nct 


vm 

'^OME,o 


- r?(T T 
9a\-^a,l -^a,2) 


(27r)3 

y)C,(r„, r„)/(m + j, o - ml, 1,0 - hl,2n„,b). 


(64) 


All of the terms in Eq. (|6^) are known, and the potential can now be calculated numerically. 
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3. TME Potentials 


We two-meson-exchange (TME) potentials consider here are the box diagrams (see Fig. ^ 
and the contact diagrams (see Fig. ||). We do not consider the crossed diagrams because 
they are not needed to restore rotational invariance, as shown in section EI3- 


4- Stretched Box Potential 


We use the Feynman rules to write the stretched-box potential shown in Fig. §(b), 


T,oif,ai 

'^TME:SB 






(Pk2m±dk 


2m 


(27r)3 J 2{2nrktA 

hm + ^-r V /(-r ^ 2 m + ^n 


(f-p hm + V (f-p hm + ^ -p V 

1/\U«/ 2M ^ “i )/I™ ^ 2/U t ay. ^aij)2n 

0{ktJ0{ktje{k+ - ktje{kt^ - ktf) 

__ 

{ktm - ktf){P- - ki^ - ky) - ml^ - - kimx) 

1 


P - /Cli - k^f -Qf -Qi 

_ Fo^Mi)^ _ 

{kti - kt,n){P~ - kii - k2m) - ml. - {ki^x - kii±)^ 

+{1 2 }. 


(65) 


To compress notation, we defined i/(= u{kif, Xif), )u = u^ku, A^), and so on. The symbol 
{1 2} means that all labels 1 are replaced with 2 and vice versa. This is a way of explicitly 

stating the potential is invariant under exchange of nucleons 1 and 2. 

We also used the following relation to simplify this expression: 


xl-2f Foif ,j Xt2, 


yt 

I AT277 


P 


O-i PK.T2i Xti f Foi f ,j Xti 


yt 

I 


T 

-L rv 


iXru 


— ■ Taf,2][Tai,l 

= [Taf,l ■ Ta^,2][Tai,l ' Ta^^2]- 


( 66 ) 


5. Mesa Potential 


We consider the Mesa potential next, Vtme:M, shown in Fig. §(a). Using the Feynman 
rules results in 


'^TME:M 


■ T„^,2][TQ,.y ■ To,_2] f (Pk2m±dk^, 


2m 


(2i) 




1 / 


7"W U ,, ,It. hm + M 


(( ^ay 2 j)lm X 2/(( Tq-j. 


2M 


raj) 


2m 


o{ktje{ktje{k+ - ktje{kt^ - k+p) 
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1 


Fa.iqff 


(Hf - Hm){P -hf- hm) - - {kif^ - 2M 

_ FaAq^Y _ 

(Hi - k2m){P~ - Ki - Km) - ml. - {ki^x_ - kii2_Y 


+{1 <—>■ 2 }. 


( 67 ) 


6 . Contact Potential 

The last potential we will calculate explicitly is the contact potential Vtme:C) shown in 
Fig. |^(a). The rules relate this potential to Vtme:M in a simple fashion. If both the initial 
and final mesons are not pions, then the potential is 


■,TOLj,ai _ 

'^TME:C ~ 


■ Taj:^2][Tai,l ' f Cpk2m±dk^.^ 


2m 


(2vr)3 J 2{27r)^ktm 

^(ktje{k+ - k+je(kt„ - fe+) 

__L 

{ktf - ktK{P- - k^j - k^m) - ml^ - {kif± - ki^±y 2M 
_ FKq^)^ _ 

{kti - k2m){P~ - kii - k^K) - ml. - {kim± - ku±)‘^ 

+{1 > 2}. 

To get the contact potential for the pions, the following change has to be made: 


[Taf,l ■ Ta.2][Tai,l ' 


r| = 3. 


( 68 ) 


(69) 


This is due the additional symmetry that the two pion vertex has. 

Before this potential is used in calculations, it must be multiplied by a factor of 2 (M+m ) 
(as discussed in section [11 E 1|) if the crossed diagrams are not included. 


7. Common Features and Simplification 

Each of the potentials in Figs. § and ^ can be written schematically in the following 
general form: 

^™E (27r)^ J 2(27r)3 i % ) 

lf{{^afM.l^aY))lrn X 2/( (Fq^ AI 2FQ, ■) ) 2m 

«/ - ml^ + bj cos(0/ - 0m)] Faf (q/)^ 

[am + bmf cos(0/ - 0m) + &m* COs(0m “ 00 ]""" 

cos(0m - 0i)l * (70) 
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We use the expression for Fa in Eq. (^ 8 ]) to get 

r (Pk2m±dkt^ 


Vtme — 


( 2 vr) 


n(j + ,+ + +N 

2(2^)3 •> > 


l/((ro^Af irc^,^)lm X 2f{{pafM.2^a^)2m 




«/ - ^af + bf COS(0/ - (j)m) O/ “ + bf COs(0/ - 0^) 


2 tiq 


"h bmf COs(0j (pm) “1“ ^mi COs(0m 0i)] 


+ bi cos{(pm - (pi) 


A^. + hi cos{(pm - (pi) 


2 tiq 


( 71 ) 


The TME potentials, like the OME potential, have exact rotational invariance about the 
z-axis, and thus conserve m. We project the potential onto states of dehnite m by setting 
(pi to zero and integrating the potential with / d(pfe^'^^f /27r. 


Now what we have to obtain the uu matrix elements. In reference shows that we can 
write 


i/((raj.AdirQA)]^^ X 2/((rQ,j.Ad2ra. 


2m 




j i^f—0 

. pif^<pn 


(72) 


Note that the C is implicitly a function of the p±, (p and A on both sides. 

Then, after a change of variables, (pm —^ (p'i, (pf — (pm (p'f, and (pf —i- (p'j + (p[, the 
potential can be written as 


_ fdk2m±k2m±dktmrn+ ,+ + +^ 

^27r)^ / ^ J 2(27r)^ / (^Im? ^2m5 0^/ : ) 


Cj^kO^af, -Ml, ^ ai]^ af, ■M.2-1^ ai) 

^277 |■27^■ d(p^ ^i(m+j)(j>'f ^i(m+k)<l>', 


/o 27r Jo 27 r 


2 / f 2 / 

ttf -rUa^ +bfCOS(pj ttf - Aaj^ + bfCOS(pf 

dm T ^m/ cos (p J T COS 

di- rria. +bi cos 0 - * - A^, + bi cos 0 - 


— 2 tiq 


2 ? 1 q 


(73) 


The azimuthal-angle integrals taken together are denoted as /, and the method for calcu¬ 
lating the (p'^ and 0^ integrals is discussed in Ref. Summarizing, we can write 


_ gHM‘^T^{A‘^ - 771^)^”“ ^ f dk2m± /^2m± dktm 

- AJ W 

^ oii]^ a; 1 -^. 21 ^ ai) 

- A^^, 6/, n/, 2n„j,, m-F j; 

dmi bmf 1 bmii (dmi 

tti - ml tti - Al bi,ni,2nai,m + k). 


fi.ktm^ktmA'fAt) 


(74) 


This potential is evaluated using the numerical integration techniques discussed in Ref. . 
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8. Further Development of the Light-front Schrodinger Equation 


The potentials derived here possess a high degree of symmetry. To solve the light-front 
Schrodinger equation efficiently, these symmetries should be explicitly exploited, as was 
done in Refs. |H, M. In addition to the invariance of the potentials under parity, there are 


additional symmetries due to the conservation of nucleon helicity and invariance under time 
reversal 
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We follow Machleidt’s approach for taking advantage of the symmetry of rotationally 
invariant potentials with helicity [^. However, since the light-front potentials derived here 
do not have full rotational invariance, Machleidt’s approach must be modified. The sym¬ 
metry properties of helicity matrix elements are rederived in Ref. |iS] without assuming full 
rotational invariance. These results allow a modified version of Machleidt’s approach to be 
combined with the exploitation of parity (using the transformation from light-front coor¬ 
dinates to equal-time coordinates) discussed in Refs. m P!P| . In particular, the potentials 
are initially calculated in the |pet, M, Ai, A 2 ) basis, although the relations in Ref. are 


used to transform to the |pet, J, M, L, S) basis to solve for the wave functions. Note that 
the potentials connect states with different J values in general. Once the symmetries have 
been explicitly expressed, we can discretize the Schrodinger equation as done in Ref. 
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Note that the transformations applied to the potential in order to simplify the solution of 
the wave functions. Once the wave functions are obtained, we may apply the inverse of the 
transformations to the wave functions to express them in the helicity basis (|p_L,p’'', Ai, A 2 )) 
or in the Bjorken and Drell spin basis (|p±,p’'",mi,m 2 )) 


III. DEUTERON BINDING ENERGIES 


The next step towards numerically calculating the bound states for these potentials is 
to choose the parameters (meson masses, coupling constants, etc.) for the potentials. We 
consider the full nuclear model where the nucleon-nucleon interaction is mediated by all 
the six mesons shown in Table |. For numerical work, use the parameters for the light- 
front nucleon-nucleon (LFNN) potential from the work of Miller and Machleidt Those 


parameters were fit for a potential that used a retarded propagator for the energy in the 
potentials. Since the potentials used in this paper have energy dependent denominators, the 
parameters must be modified somewhat. We choose to vary the coupling constant for the a 
meson. The parameters are given in table |. 

As with all the other deuteron models presented in this paper, the light-front one-meson- 
exchange (OME) potential breaks rotational invariance and causes a mass splitting of the 
deuteron states with different magnetic quantum numbers. We expect that the splitting will 
be removed somewhat by including higher order potentials. 

The first step is to determine which two-meson-exchange potentials to use. One choice is 
to use only the two-pion-exchange potentials, TPE and ncTPE, as dehned in the previous 
section. However, we expect to get better results using the two-meson-exchange diagrams 
generated by all the available mesons, including the contact diagrams for the pions, which 
we denote as the two-meson-exchange (TME) potential. In addition, we can also investigate 
the effect of leaving out the contact potentials for the pions, resulting in the non-chiral 
two-meson-exchange (ncTME) potential. 

We do not include diagrams with a contact interaction between the nucleon, a pion, 
and another meson. This is because, as mentioned in section [H E 1|, the inhnite series of 
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the box diagrams is rotationally invariant and the contact diagrams are not needed to 
achieve rotational invariance. Furthermore, they are not required to control the convergence 
of the series, since there is no strong cancellation between the contact diagram and the 
instantaneous diagrams. 

The first step in analyzing the bound states is to determine what range of fa- gives 
reasonable results. We iteratively solve for the binding energy of the deuteron, varying fa 
until the binding energy matches the physical value of the binding energy, for each of the 
potentials. The results are shown in Table We hnd that a value of fa in the range 1.2 to 
1.3 will give reasonable results for the binding energy. Note that D-state probability (about 
3%) is lower in this model than for the energy-independent light-front used in Ref. p0| , 
were a value of 4.5% is found. This is expected since the fa is greater than 1 in this model, 
meaning that the scalar interaction is strengthened relative to the tensor interaction, leading 
to a decrease in amount of the D-state present. 

We choose two values of fa, one from the low end of the range (1.22) and one from the 
high end (1.2815) for our investigations. Using two values helps ensure that our results are 
robust. 

First, we examine the bound states for fa = 1.22. The results for several different 
choices of the TME potentials are shown in Table IP and the binding energies are plotted 
in Fig. p. In addition to the two-meson-exchange potentials mentioned above, we also 
consider the vr-a plus vr-u; Mesa potential. The reason for considering this potential is 
that Carbonell, Desplanques, Karmanov, and Mathiot have shown that it helps restore 
rotational invariance of the deuteron. 

In particular, they have used manifestly covariant light-front dynamics to analyze the 
deuteron. They start with a deuteron wave function calculated in equal-time dynamics, 
then use a light-front one-pion-exchange potential (expanded to lowest order in powers of 
■^) to calculate the perturbative corrections to the wave function. They hnd that the 
resulting wave function has an unphysical dependence on the orientation of the light-front 
plane, which would manifest itself as a breaking of rotational invariance in our formalism. 
They also use the tt-ct and n-uj Mesa potentials (expanded to lowest order in powers of 
^) to calculate the correction to the wave function. When the wave function corrections 
are combined, they hnd that the directional dependence of the longest range part of the 
deuteron wave function cancels exactly. 

This implies that for our model, using the tt-ct plus tt-o; Mesa potential (which we denote 
by 7r-(cT-a;)) should partially restore the rotational invariance of the deuteron, assuming that 
the breaking of rotational invariance is due primarily to the one-pion-exchange potential. 
Note that since we solve for the deuteron wave function self-consistently and to all orders 
for our potentials, we do not expect to hnd exactly the same result as Ref. . 


The hrst thing to notice about the data in Table is that the results are essentially the 
same regardless of if arbitrary angular momentum or is used or if the potential is restricted 
to the J = 1 sector. The same result is also seen in the pion-only model . It means that 
the wave functions are numerically approximate to angular momentum eigenstates. 

Next we notice the splittings between masses and D state percentages for the m = 0 and 
m = 1 states. This implies that the states do not transform correctly under rotations. All of 
the two-meson-exchange potentials used reduce the splittings by similar amounts, by about 
60% for the binding energy and by about 70% for the percent D state. Note also that the 
mass splittings for the pion-only model were much larger pSf . 

Examining the effects of the individual two-meson-exchange potentials, we see that 7 T-{a- 
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a;)) potential does reduce the mass splitting, but it does not fully remove it. This is expected 
since the OME potential includes more than just the pion potential, and the potential is 
relativistic. 

Next, we compare the ncTME and ncTPE potentials to the TME and TPE potentials. 
The non-chiral potentials reduce the binding energy more than the chiral potentials, as we 
expected from our experience from the pion-only model. However, unlike for the pion-only 
model, we hnd that the chiral and non-chiral potential have fairly similar effects . 


Finally, notice that the mass splitting for the TPE potential is much smaller than for 
the other two-meson-exchange potentials. By itself, this does not imply that the rotational 
properties of the deuteron calculated with that potential are signihcantly better than those 
from other two-meson-exchange potentials. The individual potentials that make up the 
TME potential are fairly large in magnitude, but vary in sign. This means that using any 
subset of those potentials may result in either a larger or smaller mass splitting. In this 
case, it is smaller. To investigate this further, we examine the currents for the TME and 
TPE deuteron wave function in section |I^. 

To verify that our results are independent of the value of /o-, we recalculate the deuteron 
properties for each of the potentials with /o- = 1.2815. The results are summarized in 
Table |I^, and the binding energies are shown in Fig. The change in f^j increases the 
binding of the states, but the rest of the results are qualitatively the same. 


IV. FORM FACTORS OF THE DEUTERON 

In section 0 , we considered several different truncations of the light-front nucleon- 
nucleon (LFNN) potential and used them to obtain the deuteron wave function. In this 
section, we use those wave functions to solve for the matrix elements of the deuteron current 
operator, which is used to calculate the deuteron electromagnetic and axial form factors. 

In this section, we hrst outline the covariant theory of the electromagnetic form factors 
for spin-1 objects, like the deuteron. Then we recall the features of light-front calculations 
(including the breaking of rotational invariance) of the form factors. After that, we review 
the covariant and light-front tools for calculating axial form factors. The formalism is then 
applied to calculate the electromagnetic and axial currents and form factors for the light- 
front deuteron wave functions. 


A. Electromagnetic Form Factors 

1. Covariant Theory 

In the one-photon-exchange approximation, shown in Fig. ^ the amplitude of the scatter¬ 
ing process ed —> ed is just the contraction of the electron and deuteron currents, multiplied 
by the photon propagator. 


(75) 

where 

{p'A'\jl\pA) = Cu{p\X')^^u{p,\). (76) 
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From Lorentz covariance, parity invariance, and time reversal invariance |^, |^, 

53, we infer that the deuteron form factor can be written as 


{k',m'\J^\k,m) - 

F(fc',m')T^^e‘"(fc,m), 

(77) 

where the spin-1 polarization vectors satisfy 



e*p{k,m)e'^{k,m') = 


(78) 

0e*(fc,m)e^(fc,m) = 

m 

kfqk^ 

+ ^2 1 

(79) 

kp,e^{k,m) = 

0, 

(80) 


and the operator is given by 






(81) 


where = g^g'^ — QpQ^ is the generator of inhnitesimal Lorentz transformations, q = k' — k, 
and Fi, F 2 , and Gi are functions of the invariant mass of the photon. Some conventions 
Eg p!7| denote Gi as —F 3 . 

The Fi, F 2 , and Gi form factors are related to the deuteron charge, magnetic, and 
quadrapole form factors, denoted by Fc, Fm, and Fq respectively, by ® 00 


Fc = Fi + -rj [Fi + (1 + r])F 2 - Gi], 
Fm = Gi, 

Fq = Fi + {l + g)F2 - Gi, 


(82) 

(83) 

(84) 


where 


V = 



(85) 


Note that since both the initial and final electron and deuteron are on-shell, —> 0. We 
dehne Q‘^ = —q^. 

At zero-momentum transfer, the deuteron charge, magnetic, and quadrapole form factors 
are simply 


Fc{^) 

= 1, 

(86) 

l^A/(0) 

= = 1.71293, 

(87) 


nip 


^’e(o) 

= MlQd = 25.8525, 

(88) 


where rup is the mass of the proton, fid is the magnetic moment of the deuteron (in nuclear 
magnetons fiN = 2^)5 Qd is the electric quadrapole moment of the deuteron (in units 
of the deuteron mass, although typically it is measured in e-barns). The numerical values 


of the form factors are determined by using experimentally measured values |52, 55, 56 
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Another set of form factors are Go, Gi, and G 2 , commonly used in light-front dynamics 
r7| . They are simply related to the charge, magnetic, and quadrapole form factors. 


13 


Go = 
Gi = 

G 2 = 


Fc, 
—Fm, 
Vs p 

-Y’iFq- 


(89) 

(90) 

(91) 


The electron-deuteron cross-section is measured to determine the deuteron form factors 


for large momentum transfers. The unpolarized cross-section has the form 


da 

dfle 


CMott 


1 + 


2 Ee 

Ma 


sm 


2 


A{Q^) + 



where 


O'Mott — 


a cos Y 
2£. sin^ % 


(92) 


(93) 


and Ee is the electron beam energy, 9^ is the angle by which the electron scatters, and a is 
the hne-structure constant. The structure functions A and B can be expressed in terms of 
the charge, magnetic, and quadrapole form factors, 

2 


A = 


^2 Y 

C ' 


gVF^ + ^VFL 


B = -v{l + v)F, 


M- 


(94) 

(95) 


The three form factors Fi, F 2 , and Gi cannot be determined uniquely from the unpolar¬ 
ized scattering data since that measurement of provides only two structure functions. More 
information is needed, which can be obtained from experiments where the polarization of 
the electron and/or deuteron is measured [^, |^. The most commonly measured quantity 


is the tensor polarization observable, T 20 , 

IvFcFq + IVFq + \r]Fl^ (1 + 2(1 + 77 ) tan' 


2 


^20 — 


2{A + B tan^ ^ 


(96) 


Note that T 20 depends on the angle 6 ^- For ease of comparison, data for T 20 is usually 
presented for Of. = 70°. Alternatively, one can eliminate the angular dependence by dehning 
r 2 o, which is r 2 o with Fm set equal to zero. 


F 20 — 


InFcFq + IrfFl 


Fc + 


WF^)' 


(97) 


The extraction of the structure functions from the data is difficult. In practice, each cross- 
section measurement is performed at a different momentum transfer and angle, making it 
impossible to exactly disambiguate, for example, A and B. One needs to interpolate between 
values of B to calculate values for A, and vice versa. The Jefferson Lab t 2 o Collaboration 
used a self-consistent method for obtaining the structure functions from the scattering data 
58| . They consider a variety of theoretical models for the deuteron form factors, then fit 


the parameters of the models to the measured cross-sections. After each model is if, it is 
used as the the interpolation function to disentangle the structure functions. This two-step 
process helps minimize model-dependent effects in the values of A and B. 
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2. Light Front Calculation 


Light-front dynamics is particnlarly well snited to calculating form factors. One reason 
is that the generators of boosts in the one, two, and plus directions are kinematic, so that 
wave functions calculated with a truncated potential transform correctly under boosts. This 
feature is especially important for form factors at high momentum transfer since the wave 
functions must undergo a large boost. 

Another, more subtle, reason for using the light front is that many of the graphs which 
contribute to the current vanish identically. For example, the three lowest-order graphs for 
the current are shown in Fig. |^. The double line denotes the deuteron, and the vertex of 
the deuteron lines and the nucleon lines represents the deuteron wave function. The graph 
labeled (a) does not vanish and is calculated in section |1V A 5 . 

Fig. ^b) vanishes in light-front dynamics. To see why, we first note three facts: the plus 
component of each particle in light-front dynamics is non-negative (for massive particles 
it must be positive), the plus component of the momentum is conserved, and the plus 
momentum of the vacuum is zero. Combining these facts, we find that any vertex which has 
particles on one side and vacuum on the other must vanish. In other words, the vacuum is 
trivially empty, and no graphs couple to it. 

For Fig. 11(c), the coupling of the photon to the nucleon goes like where the 

light-front spinors are dehned in the Appendix. For /x = -(-, this matrix element is suppressed 
maximally, and thus is the “good” component of the current [^, 0, We calculate 
only J+, since it is the most stable. 

We do not consider the contribution of higher-order graphs to the deuteron current, such 
as graphs where the photon couples to a meson or a nucleon while a meson is present. 
The omission of the meson-exchange currents does not affect the rotational properties of 
the current, although it does affect the overall values of the deuteron form factors. This is 
acceptable since we are only concerned with the rotational properties in this work, not in 
the detailed results. 

The neglect of the graphs where the photon couples to a nucleon while a meson is present 
may affect the rotational properties of the current. This is because the deuteron current 
shown in Fig. ^(a) is not formally conserved |^, but current conservation is necessary 
(although not sufficient) for the current to have the correct properties under rotation. To 
construct the conserved current operator associated with a given wave function, the current 
must include diagrams that are related to the potential used to calculate the wave function. 
We expect that these diagrams are small since they contain meson propagators, and that 
neglecting them should not significantly affect the conservation of the deuteron current or 
the rotational properties of the current. 


3. Symmetries of the Electromagnetic Current 

Now, we use symmetries to relate the components of (fc',m'|J+(g)|/c,m). Although in 
light-front dynamics some generators of Lorentz transformations are dynamic, such as the 
generators of rotations from the x-y plane into the x direction are dynamic, we are free to 
boost in the plane perpendicular to the 2 ;-axis, boost in the plus direction, and rotate about 
the 2 ;-axis, since the generators of those transformations are kinematic. In addition, we will 
hnd how the states and the current operator transforms under parity and time-reversal. 

The kinematic generators allow us to choose which frame to evaluate the current in. We 
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choose the Breit frame |^, where q'^ = q~ = = 0 and q^ = Q. We also choose the plus 

momentum of the deuteron to be since Md is value of the plus momentum in deuteron’s 
rest frame. To simplify notation, we define the matrix elements of as 




(98) 


This quantity is represented with Im'm instead of Jm'm because Jm'm is used to represent the 
matrix elements of J using the instant-form spin basis. This distinction is discussed later in 
this section. 

Next, we consider rotation about the ^(-axis by tt (i?^(7r)), parity (If), and time-reversal 
(0). First note that the current operator transforms as 

= J±, (99) 

nj^n = j^, (100) 

0J±0 = (101) 

Obviously, we must use both parity and time reversal to retain the plus component of the 
current. Now look at how the deuteron states transform: 


i?2(7r)| - fc_L,m) = I-F fc_L,m)(-l)”', (102) 

0n| - = I - fc^,-m)(-l)”^. (103) 

Although the symmetry operators can introduce extra, constant phase factors, they are 
omitted here since they cancel when calculating matrix elements. These results make it 
easy to hnd what each of these symmetry operators does to the matrix element. Under 
rotations a rotation of vr about the ; 2 -axis, the current matrix element transforms to 


+T'’" 


n,,UQ) - 

Under parity followed by time reversal, the matrix element becomes 


(104) 


(105) 

(106) 


Since is Hermitian, we can also take the complex conjugate of the matrix element to get 


i(Q) 




, m 


J*(Q) 


+T'™ 


^m' V / \ 2 

The appropriate combinations of Eqs. (|104| - |107|) give []^ 




(107) 

(108) 
(109) 


The same relations also apply for matrix elements. Eqs. (|108|) and ( [109|) imply that of 
the nine possible matrix elements of J"*", there are only four independent components. We 
choose those components to be /{g, /iLi, and Iqq. It is helpful to express the matrix 
elements in a matrix to see the symmetry properties explicitly. 


-'ll -'lo -'i-i 


_ j+ j+ 

-'lO -'oo 


7i+o 


, j+ _ j+ j+ 

\ -'l-l -'lo -'ll 


( 110 ) 
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Rotational Invariance and the Angular Condition 


This is as far as we can go with light-front dynamics, bnt there should be an additional 
redundancy in our matrix elements. We have derived four independent components, whereas 
in a fully covariant framework there are only three form factors. The resolution of this 
conflict is that full rotational invariance imposes an angular condition on the light-front 
matrix elements. 

The angular condition can be found by using the deuteron polarization vectors for the 
Breit frame to calculate the current given in Eq. ( |77[) in terms of Fi, F 2 , and Gi. By 
comparing that result with Eq. ( |110|) , we obtain linear relations between the current matrix 
elements (/{^, lt_i, and Iqq) and the form factors (Fi, F 2 , and Gi) and also find that in 
general there is a deviation from the angular condition, which we denote with A, given by 

0 


A = 


-4+ + (1 + 2r/)/+ + J +1 - 2^/+, 


( 111 ) 


where rj is dehned by Eq. (pSl) . Since A vanishes when the deuteron current transforms 
correctly under rotations, we interpret A as a measure of the extent to which the current 
transforms incorrectly. 

The form factors are overdetermined by the current matrix elements, which means there 
are many different relations between them are possible. When A is zero, the angular con¬ 
dition can be used to show that all the relations are equivalent, while a non-zero A means 
that the form factors depend also on which relation is chosen. 

Since A is non-zero in general, it is important to choose the best relation to obtain the 
form factor. To do this, we classify the current matrix elements as either “good” or “bad”. 
This classification is similar to the one made for choosing which component of the current 
to use. 

In this work we are interested in the overall rotational invariance properties; the choice of 
how to relate the form factors is simply useful for comparing with other approaches. In that 
spirit, we consider four different choices [^. First, Grach and Kondratyuk (GK) derive a 
relation in which 
the equations, they obtain 


is considered bad. Using the angular condition to eliminate the for 
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Gq^GK 

Gi,gk 

G2,GK 


= 21 


(3 - 2r/)/+ + 2^1+ + /+1 
1 


11 




no ) 




Brodsky and Hiller (BH) use a prescription where is bad, which results in 

1 


Gq^BH — 
Gi^bh = 
G2,BH = 


3(1 + 277) 
2 


(3 - 2 ,,)/+ + Sd^il+ + 2 ( 21 , - 1 ) 7+1 


1+277 
2 V 2 

3(1 + 277) 


7„+„ - 7i+_. + (2i, - 1) 


It 


v+At - ’At - (1 + 


( 112 ) 

( 113 ) 

( 114 ) 

( 116 ) 

( 116 ) 

( 117 ) 
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Frankfurt, Frederico, and Strikman (FFS) start by analyzing the deuteron current using 
an equal-time spin basis, where the current matrix elements are <7^/^ [^. Their prescription 
uses the Cartesian basis to hnd that is the bad current. Upon transformation to the 


spherical basis and use of the Melosh transformation [^, to relate 
hnd 


to 


they 


Go,FFS 

Gi,ffs 

G2,FFS 


3(1+ h) 

1 


{2r] + -h {2r] + 


l + r] 

V2 


3(1+ h) 


/+ I /+ _ /+ _ 

“ “ FSj 

Vlu + SyU-fJi - -(v + 2)y_i 


(118) 

(119) 

( 120 ) 


Chung, Polyzou, Coester, and Keister (CCKP) choose the canonical expressions for the 
form factors Gq, Gi, and G 2 in terms of the equal-time current as the starting point 0,0- 
They use rotations and the Melosh transformation to express the equal-time current matrix 
elements in terms of the light-front current matrix elements, and hnd that 


Go, CCKP 

Gi,cckp 

G2,CCKP 


1 

6(1+ r/) 

Gi,ffs, 

G2,FFS- 


(3 2r7)(/j^-|-/ q())-|--|-( 4?7 l)7{*Li 


( 121 ) 

( 122 ) 

(123) 


5. Impulse Approximation on the Light Front 


We now are ready to relate the deuteron wave function to the current expressed in 
Eq. (pSl). The deuteron wave functions solved in section II are used. We use the repre¬ 
sentation of the wave function in the spin \p±,p'^, 1711 , 1712 ) basis, as discussed in Sec. fimj . 


Note that these spins are expressed in the usual Bjorken and Drell representation . By 


inserting a complete set of states into Eq. 
explicit, we get 


and making the momentum of particle 1 






2(p+,px + 


07 

2 ’ 




J{S)i07) 


+ 07 

P AP7 - 


07 


P ,P7- —,1711,7712 


07 

-, m , , 

2 ' 


(124) 


where the spin directions of particles 1 and 2 are labeled irti and m 2 , respectively. Also, 
since the deuteron is an isoscalar combination of nucleons, the isovector component of the 
nucleon current does not contribute and the isoscalar nucleon current is the same for both 
nucleons. This allows us to simply double the isoscalar current of particle 1 instead of using 
the isoscalar currents of both particle 1 and 2. 

We want to boost the deuteron wave functions to the rest frame, since that is where the 
wave functions are calculated. Note that the boosts in the x-y plane transforms a general 
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vector (fc^, /c+) by fc_i_ ^ fc_i_ — q±^, leaving k~^ unchanged. This means that the boost that 
puts the deuteron in its rest frame transforms the individual nucleon momentum by 


P±- Y 


(125) 


where we use the Bjorken a:-variable x = ^. Note that the spin labels do not change. This 
is because the boost is in the perpendicular direction, so the spin labels (dehned to point in 
the z direction) are not affected. 

This allows Eq. (|124|) to be written as 






2(p+,P± + 




+ 

p\px - 


- (1 - a:)^,mi,m2 rri^ . 


(126) 


We have dropped explicit mention of the deuteron’s momentum. 

Before we continue, we note that we can write the nucleon current matrix elements using 
the u spinors. 


P^,P± + 


'^(5)(9±) 


+ 

p\p± - Y^^^ 


= Ubd {p^,P± + ) T 5 )(q±)MBD (p^,P± “ ) , 


'(5) 


q± 


(127) 


where we have used the “BD” subscript to denote these as Bjorken and Drell spinors. 
Since the nucleons are on-shell, as they must be in a Hamiltonian theory, we interpret 

wbd (p+,p±,m) as mbd (px = P±,x,Py = P±,y,Pz = where p" = ^ . 

However, the matrix elements can be calculated easily if we convert the Bjorken and Drell 
spinors to light-front spinors. This transformation is formally accomplished with the Melosh 
transformation Rm, given by |^, p!7| , 

M — icr ■ {n X p±) 


Rm — 


p' 


( 128 ) 


\/(p+ + M)2 + pi 

where n points in the direction of the light front. Note that the Melosh transformation is 
unitary. 

This transformation converts a Bjorken and Drell spinor to a light-front spinor in the 
following manner, 

ULF{p"",P±,m) = J2idiM)m,m'UBD{p'^,P±,m'). (129) 


Using this transformation, we rewrite Eq. 


as 


= J d^Pxdp-^ ^ p-^,p± + {l-x)Y^^l^m2'^{RM)m';,m[ 


2 u 


LF 


L^i 

q± 


P ,P± + —) ^( 5 )( 9 ±)wlf (p - —,mi 


q± 


IZCfilljm',"..,,, - (1 ™) 


( 130 ) 

























6 . The Nucleon Form Factors 


From Lorentz covariance, parity invariance, and time reversal invariance |^, 
isoscalar part of the nncleon cnrrent can be expressed as 

= YFl^''(g) + i^FP(q). (131) 

Taking the matrix elements of with the light-front spinors gives 
u(k', A')J^j(g±)u(k, A) = Fl^\g)uLF(k', A')'y~^ULF(k, A) 


51 , the 


+FF(<!)^ULF(k', \')7*7[uLF(k, A), 


(132) 


where we have parameterized 


= k+, 

k± = p± + —, 


k'± = P±- 




(133) 

(134) 

(135) 


We nse the representation of the spinors given in the Appendix to simplify the spinor matrix 
elements appearing in Eq. |132|. First consider 


UFF{k',A')-f^UFF{k,A) = 


Mk+ 


Xlf,a' 


X 


\m - ■ k'^)A+ + k+A_ 

A_{M + a"*" ■ fc"*") -|- A+/c 


7 

Xlf,a, 




Next, 


mlf(^', A07^7±wlf(/s, A) = 


Mk+ 


Xlf,a' 


X 


\m - ■ fc'^)A+ + fc+A_ 

A_{M + ■ k^) + A+k+ 


7 7 
Xlf,a, 


k+ 


A 


— ^Xlf,a'7 7±Xlf,a, 


Combining these, we hnd M, 12 


{k',\'\Jt7,M\k,\) = ^aI [fkh) - 7) 

We can rewrite Eq. (|141|) as 

(4',AV(t)(<l±)|i,A) = -F7(<;)(i',AVas)(<;i)|i.A> + 

FT{q){k',\'\J*,^(qj.)\k,\), 




X\ 


(136) 

(137) 


(138) 

(139) 

(140) 

(141) 


( 142 ) 
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( 143 ) 


which, when inserted into Eq. (|130|) , gives 






(2)m',m \ 


where we dehne 

I 


{i)m' ,m 


(q) = / d^p^dp-^ Y. 


(m'|p+,p_L + (1 - x)^,m[,m2) 


,mi ,7712 

q± 

2 


2(p+,px + 


'{iS) 




{P\P± - (1 - x) — ,mi,m2\m), 


(144) 


for z = 1,2. Note that both i^ust satisfy the same equations 

as rniQ) does. In particular, this means that the angular condition should apply to 
'^( 2 )m'm('?) Independently, so we consider the deviation from the angular 
condition for each. 

There are many parameterizations of the isoscalar nucleon form factors Fi^\q) and 
F 2 (q)- Since the measurement of the electron-nucleon cross section is difficult, the data 
have large errors and are consistent with several different models of the nucleon form fac¬ 
tors. Some of the models representative of those proposed in the literature are: the dipole 
model, £t 8.2 of Hohler |Q, Gari:1985 model 3 of Gari:1992 |^, best £t for the 
multiplicative parameterization of Merged [^, and model DR-GK(l) of Lomon 
Fi^\q) and F^i^) form factors for each of these models are shown in Fig. |T0. 


69 . The 


We can relate the isovector form factors to Gep, Gmp, Geui and Gmu-, the proton electric, 
proton magnetic, neutron electric, and neutron magnetic form factors, respectively, with 


70 


Ff) = 
= 


Gep + Geu + T {Gmp + Gmu) 

2(1+ r) 

—Gep — Geu + Gmp 


G 


Mti 


( 145 ) 

2(1+ r) ’ 

where r = The value of r is approximately 1 at a momentum transfer of about 5 GeV^, 
the upper range of momentum transfers that we consider. Since the overall magnitudes of 
the form factors are similar at this momentum transfer, it is important to measure each of 
the form factors with the same accuracy and cover the same range of momentum transfers. 
Gurrently, the most poorly known form factor is Geu, both in terms of the magnitude of 
the error and in the number of data points 


69 


In section 0, we will hnd that for momentum transfers greater than about 2 GeV^, the 
spread in the values of the deuteron form factors due to the breaking of rotational invariance 
on the light front is smaller than the spread in values due to using the various nucleon form 
factors. It is uncertainty of the nucleon form factors, not the use of the light front, that 
limits the accuracy of the deuteron form factors at large momentum transfers. Only more 
accurate measurements of the nucleon form factors, especially Geu, will allow for more 
accurate deuteron form factor calculations. 


B. Axial Form Factors 

The formalism used for the axial current and form factor is very similar to that used for 
the electromagnetic current and form factor. Thus, most of the discussion from the previous 
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section carries over here with only slight modihcations. We highlight only the differences. 

The derivation of the symmetries of the axial current matrix elements is almost the same 
as in section |1V A 2| , with the exception that under parity, the axial current picks up a 
negative sign 


nyn = -j^. 


( 147 ) 


When this is propagated through the algebra, we find that 


( 148 ) 

( 149 ) 


Eqs. ( |148|) and (|149|) imply that of the nine possible matrix elements of , there are only 
two independent components. We choose those components to be .f( 5 )ii, and If is 

helpful to express the matrix elements in a matrix to see the symmetry properties explicitly. 




We have derived two independent components, but an analysis of the covariant theory 
shows that only one deuteron form factor (T^) contributes for the plus component of the 

0 , 


/ 

J+ 

H5)ll 

T+ 

H5)10 




_ J+ 

H5)10 

0 


(150) 

V 

0 

J+ 

H5)10 




axial current 


This implies that the requirement of full rotational invariance imposes an 


angular condition on the light-front axial current matrix elements. The deviation from the 
angular condition, denoted by A, given by |TT[| , 


^ _ V^ j+ _ j+ 

^ “ 2 (5)11 -'(5)10- 


(151) 


Since the deuteron axial form factor is overdetermined by the current matrix elements, 
we need to classify the current matrix elements as either “good” or “bad” to eliminate 
ambiguity. We consider two such choices. 

Frankfurt, Frederico, and Strikman (FFS) find that the is the bad matrix element 


T^ . After transforming to the spherical basis and using the Melosh transformation, they 


find that 


Fa = 


2(1 + r,) 


+ 


(5)10 


( 152 ) 


Frederico, Henley, and Miller (FHM) use the behavior of the matrix elements in the non- 
relativistic limit to determine that the bad element is 


111. This means that 




(153) 


The current matrix elements are calculated using the nucleon axial current. The general 
form of nucleon axial current is given by ||11||, 


^5).. = 


n 

P- 


( 154 ) 
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Since we choose to /i = + and work in the Breit frame, where = 0, Eq. (|154|) reduces to 

Jhn = ( 188 ) 

can be computed using expressions given 


The light-front spinor matrix elements of Eq. 
in the Appendix. The result is 




(fc',AV(t)n(9±)|fc,A) = 


M 


(156) 


We consider only one model for the nucleon axial form factor since the deuteron axial 
current has such a simple dependence on it. We choose to use the dipole model 


Fa{Q^) = 


Fa{0) 


1 + 


m 


2 5 


(157) 


where Ma is the axial mass. For our calculations, we use the value for the axial mass 


determined by Liesenfeld et al. [71 


V. RESULTS FOR THE FORM FACTORS 


We use the deuteron wave functions obtained for the light-front nucleon-nucleon potential 
in section II to calculate the deuteron currents and form factors. This gives a solution where 
light-front dynamics is used consistently throughout. For the potential, we choose the light- 
front nucleon-nucleon potential with = 1.2815. We have verihed that the results do not 
change signihcantly when = 1.22 is used. 

Figure IT show the currents and the associated angular condition for given by 
Eq. (|144|) , for several different deuteron wave functions. Results are shown for the wave 
function from the OME, OME-I-TME, and OME-I-TPE potentials (calculated in section II), 
and the parameterization of the deuteron wave function for the energy independent Bonn 

HI' 


potential 


The currents matrix elements (but not A) are approximately the same re¬ 


gardless of which wave function is used. This consistency is important, since it verihes that 
the gross features of all the models are the same. 

We hnd that A for is much smaller than the largest matrix elements when using the 
OME wave function. This means that the current transforms very well under rotations. 
This is somewhat surprising, since we found earlier that the binding energies for the OME 
wave functions have a large splitting, indicating that OME wave functions transforms poorly 
under rotations. 

Comparing the current calculated with the OME wave function to those calculated with 
other potentials, we hnd that for momentum transfers of more than 1 GeV^ the OME 
current has the best transformation properties under rotation of all the /A currents shown. 

For smaller momentum transfers, the transformation properties of the Bonn and 
OME-I-TME wave functions are the best. This is expected, since in the limit of no momen¬ 
tum transfer, the current is simply the overlap of deuteron wave functions, (m'jm). 

If the initial and hnal states have the same mass, the matrix element is simply 5m'm, which 
satishes the angular condition. However, if the states do not have the same mass (which 
implies that m' ^ m), there will be a non-zero overlap between the two states, which violates 
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the angular condition. Since the masses of the deuteron states are exactly the same for the 
Bonn wave function, and approximately the same for the OME+TPE wave function, they 
have a small A at low momentum transfer. However, the OME wave functions, having the 
largest mass splitting, have the largest A at low momentum transfers. 

Figures and show the current matrix elements and the angular condition for 
and respectively. The general features of these figures are the same as in Figure |^, 
with one important exception. In both hgures, the A for the OME wave function has about 
the same magnitude as the A’s for the other wave functions. This means that the rotational 
properties of and 


currents are approximately the same regardless of which wave 


function is used. This is a surprising result since it indicates that the rotational properties 
of the current depend more on how the current is constructed than on which wave function 
is used. 

In Fig. 1^, we hnd that the magnitude of A is almost the same as the magnitude of 
the largest matrix element of This means there is a large deviation from the angu¬ 
lar condition, and that form factors calculated with this current may depend strongly on 
which matrix element is chosen as “bad”. We show below that this is not the case for the 
electromagnetic form factors. 

We hnd that A is much smaller than the largest matrix element of the axial currents 
shown in Fig. for most values of momentum transfer. This means that the deuteron axial 
form factor will be essentially independent of which matrix element is chosen as “bad”, 
except for within the range of 1.5 to 2 GeV^. 


Now we combine the two parts of the electromagnetic current, and 


/( 2 ), with the 


nucleon form factors Fi and F 2 to get the total current. Figure |T^ shows the currents for 
and as well as the sum, /■•■. The Gari;1985 nucleon form factors are used 

|65| . We hnd that Fil^-^ gives the largest contribution to the total current, and because A is 
small for A is also small for the total current, meaning that the total current transforms 
well under rotations. Thus, in spite of the fact that A is approximately the same size as 
the current matrix elements for the deuteron electromagnetic form factors should not 
depend too strongly on the choice of the “bad” matrix element. This is especially true for 
the form factors calculated with the OME wave function. 

We calculate the form factors A, B, T 20 , and Fa using the OME wave function, and 
show the results in Fig. [T^. In general, the form factors do not depend strongly on which 
matrix element is chosen as “bad”, in agreement what what we predicted in the previous 
paragraph. The dehnitions of the “bad” matrix elements are given in sections [1VA4| and 
[tVB| . The only exception is for the B form factor, and to a lesser extent the Fa form factor, 
near where they crosses zero. This is not too surprising, since a small constant shift in any 
function near a zero-crossing has a large effect in a logarithmic plot. Also, we note that the 
FES and GGKP choices of the “bad” matrix element give the same value for B. This is a 
consequence of Eqs. (|122D and ( |123| ). 

We also use the OME-I-TME wave function to calculate the form factors A, B, T 20 , and 
Fa, which we show in Fig. |^. We argued earlier that the these electromagnetic form factors 
depend more strongly on which matrix element is chosen as “bad” that those calculated with 
the OME wave function, and that dependence is clear in this hgure. At low momentum 
transfers, the dependence on the change is fairly small, but as the momentum transfer 
increases, so does the dependence. The axial form factor is not affected as strongly, primarily 
because each wave function generates an axial current which violates the angular condition 
by approximately the same amount. 
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Since there are many different models of the nucleon electromagnetic form factors, we 
calculate the deuteron electromagnetic form factors using each of them to see what effect 
the differences have. The results are shown in Fig. Q. At low momentum transfers, all the 
nucleon form factors give close to the same results. However, when the momentum transfers 
is large, we find a large spread in the values due to nucleon form factors. In fact, this spread 
is larger than the spread of values obtained from using different “bad” matrix elements with 
the OME+TME wave functions. In other words, in order to obtain accurate results for 
momentum transfers over 2 GeV^, it is more important to determine which nucleon form 
factor to use than when “bad” matrix to use. 

Finally, in Fig. |T^, we compare the A, B, T 20 , and Fa form factors for the OME and 
OME+TME wave functions to experimental data. The “bad” component was chosen ac¬ 
cording to FFS, and the nucleon form factors of Lomon were used for A, B, and T 20 , while 
the Liesenfeld axial nucleon form factor was used for Fa- The data for A is from: Buchanan 

Abbott et al. 


et al. 72 , Elias et al. [731, Galster et al. 74 , Platchkov et al. [75 


76 , and 


72 


77 


Alexa et al. [|I| ; the data for B is from: Buchanan et al 
et al. 1^; and the data for T 20 is from: Schulze et al. 

|8T| , Garcon et al. , Ferro-Luzzi et al. [^ , Bouwhuis et al. [Q , and Abbott et al. [^ 


79 


, Auffret et al. 
Gilman et al. 




, and Bosted 
Boden et al. 


There is a rather large difference between the form factors calculated with the OME and 
OME+TME wave functions. This difference is due primarily to the fact that the OME 
wave functions are more deeply bound than the OME+TME wave functions, and it can be 
reduced by choosing a different sigma coupling constant f„ for the OME and OME+TME 
potentials. However, for our analysis of rotational invariance, it is important to keep 
fixed. 

The difference between the calculated form factors and the data is also quite large. This is 
not unexpected, since in our model of the current, meson exchange currents are not included. 
It is known that these have a large effect on the form factors at large momentum transfers 
0 P6| , ^ Including these effects would bring the form factors into better agreement with 


the data. However, we emphasize again that agreement with the data is not a priority of this 
work. Our goal is to gain a better understanding of the breaking of rotational invariance by 
the light front, and how to restore that invariance. Only after we have that understanding 
can we pursue accurate calculation of the form factors with light-front dynamics. 


VI. CONCLUSIONS 


The issue of rotational invariance in light-front dynamics must be addressed before one 
attempt to use light-front dynamics for high-precision calculations. In this paper, we have 
sought to find ways to quantify the level to which rotational invariance is broken. In addition, 
we have investigated potentials of different orders, since we know that if all orders of the 
potential were included, rotational invariance would be restored. We have used light-front 
dynamics to obtain new light-front nucleon-nucleon one-meson-exchange (OME) and two- 
meson-exchange (TME) potentials. 

In section H, we derive OME and TME potentials for a model Lagrangian for nuclear 
physics which includes chiral symmetry. These potentials are then used in section |T1| to 
calculate the binding energy and wave function for the m = 0 and m = 1 states of the 
deuteron. We find that the splitting between the m = 0 and m = 1 states was smaller for 
the OME+TME potential as compared to the OME potential gives. 

In section |^, the wave functions obtained in section H are used to calculate the form 


34 



















































































factors of the deuteron using only light-front dynamics throughout. In light-front dynamics, 
there are four independent components of the deuteron current. However, the requirement of 
rotational invariance introduces an angular condition that the four components must satisfy, 
reducing the number of physically independent components to three. The deviation of the 
calculated current components from the angular condition is denoted by A. We hnd that A 
is very small for the deuteron wave functions calculated with the OME potential. This is an 
important result, since it means that although in principle the light-front calculation of the 
deuteron current does not transform correctly under rotations, in practice it does quite well. 
The smallness of A means that any reasonable prescription for eliminating the dependent 
component of the current gives essentially the same results; the uncertainty introduced by 
the various nucleon form factors is much greater. 

We also found that A is signihcantly larger when the TME potentials are used. Since 
the results in Refs. O, O, ^ indicate that the rotational properties of the TME wave 


function are better that for the OME wave function, we interpret the increase in A as an 
indication that extra diagrams need to be included in the current calculation to restore 
rotational invariance. That is, having a wave function with good rotational properties is not 
sufficient to obtain matrix elements of the current operator with good rotational properties. 
To verify this, the contribution to the deuteron form factors due to the photon coupling to 
an intermediate two-meson-exchange type of diagram must be calculated. This may restore 
the rotational invariance of the deuteron form factor when the wave function is calculated 
using the two-meson-exchange potentials. 
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APPENDIX A: NOTATION, CONVENTIONS, AND USEFUL RELATIONS 

For a general four-vector a with components [aP, a}, a?^ a^) in the equal-time basis, we 
dehne the light-front variables 

= a° ± aP, (Al) 

a_L = (a\a^), (A2) 

so the 4-vector can be denoted in the light-front basis as 

a = (a"*", a“, a_i_). (A3) 

Using this, we hnd that the scalar product is 

a ■ b = = - {a^b~ + a~b'^^ — • 6j_. (A4) 

This dehnes with = g_^ = 1/2, gn = g 22 = —1, and all other elements of g 
vanish. The elements of g^'^ are obtained from the condition that g^'^ is the inverse of g^y, so 
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9 °‘^gi 3 \ = Its elements are the same as those of except that g ^ =2. Thus, 

= 2a^. (A5) 

and the partial derivatives are similarly given by 

= 2d^ = 2^. (A6) 

^ dx^ ^ ^ 

To hnd the physical consequences of this coordinate system, consider the commutation 
relations \p^,x^] = ig^'^, which yield 

[p^,x^] = 2i, (A7) 

[P±^^±] = (A8) 

with the other commutators equal to zero. This means that x\ is canonically conjugate 
to p\, and is conjugate to Since x~^ plays the role of time (the light-front time) 
in light-front dynamics, and p~ is canonically conjugate to x"*", this means that p~ is the 
light-front energy and that the light-front Hamiltonian is given by P~. 

In any Hamiltonian theory, particles have the an energy dehned by the on-shell constraint 
= m?. This implies that the light-front energy of a particle is 

I u2 


The independent components of the 
denoted by 


For dealing with spin, we require 


The Bjorken and Drell convention for the gamma matrices is used in this work. They 
specify that 


k- = 


III ~\~ ^J_ 

k+ ' 


(A9) 


momentum can be written as a light-front three-vector 


^LF = {k'^,ki_). 

the Pauli sigma matrices, which are 

(Ciprpy-p) 


(AlO) 


(All) 


k = P 

7 = i3ol 

7 ® = 



The spin matrices S'* then are 

E* 



(A12) 

(A13) 

(AM) 

(A16) 

(A16) 
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Using S, we can express the helicity operator as H = p ■ S, which has eigenvalues ±1. This 
is useful since the helicity is invariant under rotations. 

It is useful to define the spinor projection operators A± by 


A± = ^7^7^ = ^7V = ± «^)- 

These satisfy the requirements for projection operators, 


(A17) 


A+ + A_ = 1, (A18) 

(A±)2 = A±, (A19) 

A±A:p = 0. (A20) 

We summarize the effect these projection operators have on the gamma matrices: 


A±7° = 7°A:^, 

A±7=^ = 0 = 7=^A:^, 
A±7^ = 7=f = y^Azp, 
A±7'^ = 7'^A±, 


and under conjugation, 


= A:,. 


(A21) 

(A22) 

(A23) 

(A24) 


(A25) 


The spinors used by Bjorken and Drell polarized in the z direction, and so the 

x’s are chosen to be 


X+l/2 — 
X-1/2 = 

The light-front spinors are dehned to be 0 

1 


1 

0 

0 

1 


(A26) 

(A27) 


ui^p{k, X) = 


Xlf,x = 


VMl^ L' 
1 

xfAW 

- f Xx\ 


MA_ -|- {k~^ + Q:"'' ■ fe^)A+ 
A_(M -|- cx^ ■ fc-*-) -I- A+fc’*' 


0 


Xlf,a 

(A28) 

Xlf,a, 

(A29) 


(A30) 


where xa is the usual Pauli spinor, and the A± are the spinor projection operators dehned 
in Eq. (|A17|) . We hnd that 


upF{k,X) = 


x/mB 

1 

x/MB 


Xlf,a 

Xlf,a 


A+M + A_{k++ ■ k^) 

\m - cx^ ■ k^)A+ + fc+A_ 


(A31) 

(A32) 
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Note that these spinors are normalized to satisfy uijpik, X')uijp{k, A) = 5\i\. 

For helicity spinors, we choose the eigenvectors of the helicity operator (cr ■ p) as the y’s. 
In particular, {p ■ T:)u{p, A) = hu{p, A), where h = 2A. This choice allows us to write 


and 


u{p, A) 


W 

2M[f 


x\(p), 


Xa(p) 


1 g^g+10/2 j 


a h = +1 
if h = -1 


(A33) 


(A34) 


where C 2 = cos |, S 2 = sin f, / = ^, and h = 2A. 

When there are two fermions in the center-of-momentum frame, we can dehne 0 = 0i 
and 6 = 6 i and for particle two 02 = + 0 and 62 = — 6. This means that 


u{pi, Xi) 


W / 1 
2M [fi 


Xi.y(p), 


where i = 1, 2 and 


Xi,Ai(p) = XAi(p), 

X2,A2(p) = W-A2(p)- 


(A35) 


(A36) 
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TABLE I: The parameters for the tt, 77 , p, uj, 6, and a mesons. For the meson type, “iv” and 
“is” stand for isovector and isoscalar, while “ps”, “v”, and “s” stand for pseudoscalar, vector, and 
scalar. The a meson is also known as /o(400 — 1200), and the 6 meson is the ao(980). 


Meson 

type 

mass [MeV] 

A [GeV] 

n 

477 

/ 

.9 

vr 

iv,ps 

138.04 

1.2 

1 

14.0 


V 

is,ps 

547.5 

1.5 

1 

3.0 


P 

iv,v 

769. 

1.85 

2 

0.9 

6.1 

(jj 

is,v 

782. 

1.85 

2 

24.5 

0.0 

5 

iv,s 

983. 

2.0 

1 

2.0723 


a 

is,s 

550. 

2.0 

1 

faX 8.9602 



TABLE II: The values of fa required to give the physical value of the deuteron mass for a given 
potential and state with = m. The percent of the wave function in the D-state and in the J = 1 
state are also shown. 


Potential 


u 


% D state 

%J = 

1 


m =0 

m=l 

Diff 

m =0 

m=l 

m =0 

m=l 

OME 

1.2407 

1.2125 

0.0282 

2.87 

3.55 

99.99 

99.97 

OME 

1.2829 

1.2819 

0.0010 

2.96 

3.23 

99.99 

99.97 

+TPE 

OME 

1.2968 

1.3079 

- 0.0111 

2.95 

3.28 

99.99 

99.96 

+TME 

OME 

1.3064 

1.3121 

-0.0057 

2.99 

3.16 

99.98 

99.97 

+ncTPE 

OME 

1.3198 

1.3397 

-0.0199 

2.98 

3.21 

99.98 

99.96 

+ncTME 
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TABLE III: The values of the binding energy, percentage of the wave function in the D state, and 
the percentage of the wave function in the J = 1 state for the m = 0 and m = 1 states for different 
potentials. The a coupling constant factor is = 1.22. 


Potential 

Binding Energy (MeV) 

D state (%) 

J = 

1 (%) 


m=0 

m=l 

Diff 

m=0 

m=l 

m=0 

m=l 

OME only 

-1.7653 

-2.4200 

0.6547 

2.73 

3.61 

99.99 

99.96 

OME 

+7r-(cr-w) Mesa 

-1.9236 

-1.7021 

-0.2215 

2.80 

3.38 

99.99 

99.96 

OME 

+ncTME 

-0.4948 

-0.2646 

-0.2302 

1.97 

1.76 

99.99 

99.98 

OME 

+ncTPE 

-0.6620 

-0.4825 

-0.1795 

2.16 

2.09 

99.99 

99.98 

OME 

+TME 

-0.7861 

-0.6060 

-0.1801 

2.25 

2.31 

99.99 

99.97 

OME 

+TPE 

-0.9981 

-0.9155 

-0.0826 

2.42 

2.57 

99.99 

99.98 


TABLE IV: The values of the binding energy, percentage of the wave function in the D state, and 
the percentage of the wave function in the J = 1 state for the m = 0 and m = 1 states for different 
potentials. The a coupling constant factor is fa = 1.2815, distinguishes this table from Table B 


Potential 

Binding Energy (MeV) 

D state (%) 

J = 

1 (%) 


m=0 

m=l 

Diff 

m=0 

m=l 

m=0 

m=l 

OME only 

-3.3500 

-4.4546 

1.1046 

3.09 

3.97 

99.99 

99.96 

OME 

+'K-{a-uj) Mesa 

-3.6331 

-3.2408 

-0.3923 

3.10 

3.85 

99.99 

99.95 

OME 

-|-ncTME 

-1.3766 

-0.9901 

-0.3865 

2.67 

2.64 

99.99 

99.97 

OME 

-|-ncTPE 

-1.6532 

-1.4693 

-0.1839 

2.81 

2.88 

99.99 

99.97 

OME 

-FTME 

-1.8617 

-1.6032 

-0.2585 

2.85 

3.05 

99.99 

99.96 

OME 

+TFE 

-2.1915 

-2.2137 

0.0222 

2.95 

3.23 

99.99 

99.97 
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FIG. 1; The first several terms of the full kernel for the Bethe-Salpeter equation of the nuclear 
model with chiral symmetry. 



{a) {b) (c) 


FIG. 2; The non-vanishing diagrams for pion-nucleon scattering at threshold: (a) (b) 

and (c) Me- The mesons here are pions. 



FIG. 3; The two diagrams which contribute to the OME potential for each meson. 
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FIG. 4: The TME potentials for (a) VtmE:M (the Mesa potential), (b) TtmE;SB (the stretched box 
potential), (c) TtmE:SBI (the stretched instantaneous potential), and (d) VtmE:SBII (the stretched 
double instantaneous potential). Note that the graphs on the right side are obtained from the 
graphs on the left side by 1 2. 
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FIG. 5: The TME potentials that include the contact interaction for (a) VtmE:C (the contact 
potential), (b) FrME:SBC (the stretched contact potential), (c) VtmE;SBIC (the stretched instanta¬ 
neous contact potential), and (d) VtmE;SBCC (the stretched double contact potential). Note that 
the graphs on the right side are obtained from the graphs on the left side by 1 2. 
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Binding Energy [MeV] 



Potential Used 


FIG. 6; The values of the binding energy for the m = 0 and m = 1 states for different nucleon- 
nucleon light-front potentials. The a coupling constant factor is = 1.22. 




Potential Used 


FIG. 7; The values of the binding energy for the m = 0 and m = 1 states for different nucleon- 
nucleon light-front potentials. The a coupling constant factor is = 1.2815. 



FIG. 8; The Feynman diagram for one-photon-exchange electron-deuteron scattering. 
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and 


FIG. 9: The lowest-order graphs which contribute to the deuteron current matrix element. 



FIG. 10; The Fi and —F 2 isoscalar nucleon form factors for six different models: the dipole model, 
Hohler, Gari:1985, Gari:1992, Merged, and Lomon. 






FIG. 11; The matrix elements of the component of the electromagnetic current which 

multiplies the nucleon Fi form factor, calculated with the wave function from the (a) OME, (b) 
OME+TME, (c) OME+TPE, and (d) Bonn potentials. 
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FIG. 12; The matrix elements of the component of the electromagnetic current which 

multiplies the nucleon F 2 form factor, calculated with the wave function from the (a) OME, (b) 
OME+TME, (c) OME+TPE, and (d) Bonn potentials. 
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FIG. 13: The matrix elements of the deuteron axial current including the nucleon axial 

form factor, calculated with the wave function from the (a) OME, (b) OME+TME, (c) OME+TPE, 
and (d) Bonn potentials. 
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FIG. 14: The matrix elements for ^m'm calculated with the OME wave 

functions. The Gari;1985 nucleon isoscalar form factors are used for Fi and Fo. 
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FIG. 15; The form factors A, B, T 20 , and Fa calculated using the various choices of the “bad” 
matrix element. The OME wave function is used, along with the Lomon nucleon form factors for 
the electromagnetic form factors, and the Liesenfeld nucleon form factor for the axial form factor. 
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FIG. 16: The form factors A, B, T 20 , and Fa calculated using the various choices of the “bad' 
matrix element 


IV B 


The dehnitions of the “bad” matrix elements are given in sections IV A 4 and 
The OME+TME wave function is used, along with the Lomon nucleon form factors for the 


electromagnetic form factors, and the Liesenfeld nucleon form factor for the axial form factor. 
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FIG. 17; The electromagnetic form factors A, B, and T 20 calculated using the various choices of 
the nucleon isoscalar form factors. The OME wave function is used, along with the FFS choice of 
the “bad” deuteron current matrix. The axial form factor is not shown since its dependence on 
different form factors is trivial. 
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FIG. 18; The A, B, T 20 , and Fa form factors for the OME and OME+TME potentials, along with 
data. See the accompanying text for an explanation of the data. 
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